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Abstract 

We derive the QCD sum rules for the vector and scalar meson mixing in nuclear medium, 
using a two quark interpolating field for both mesons. Modeling the mixing via a nucleon hole 
contribution with known coupling constant, the sum rule can be used to determine the overlap of 
the interpolating field and the scalar meson. In the 1=0 channel, we find a stable Borel curve and 
an overlap that is about 10% of the corresponding value in the pseudo scalar or vector channel. 
The sum rule in the 1=1 channel is less reliable but also consistent with a small value for the 
overlap. These results suggest that both the a and ao have a small two quark and thus probably 
a large tetraquark components. We discuss the possibility of observing these scalar mesons from 
vector mesons emanating from the nuclear medium. 
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I. INTRODUCTION 



The properties of vector meson immersed in nuclear medium is expected to change due 
to partial chiral symmetry restoration and/or many body nuclear effects. Expected changes 
ranee from the decrease in the mass 111, M, increase in thewidthgHSHg, and appearance 
of new structures [1, 0] , every aspect of which if measured will lead to a new understanding of 
strong interaction. Several experiments are reporting preliminary results 10, 11, 13], which 
indeed hints to nontrivial changes and new structures appearing at nuclear matter. With 
further refined experiments and with exclusive measurement of final state particles, the 
background could be substantially reduced and vector meson kinematics controlled. Then 
detailed information on the changes of vector meson properties inside nuclear matter can 
be extracted, which will then serves as a basis for understanding; symmetry restoration in 
QCD, generation of hadron masses, and QCD phase transition [13j. While the prospects of 
looking directly into the nuclear matter through vector meson seems so exciting, a careful 
model independent theoretical analysis on the vector meson properties at nuclear matter 
have to be carried out simultaneously before any conclusion on medium effects can be made. 

An important feature of the vector meson when immersed into the nuclear matter is its 
coupling to the scalar meson with the same isospin. Therefore, the longitudinal mode of uj 
meson will couple to a, and that of p to a . The coupling is possible because the nuclear 
matter expectation value of operators with Lorentz index (spin) can be non vanishing. 
Therefore, in isospin symmetric nuclear matter, mesons with different spins can couple 
as long as their parity and isospin are identical. The coupling between the scalar and 
the vector meson have been extensively studied in the Walecka model in the isospin 
channel 14, 15| . and in the isospin 1 channel 16, 17]. On the other hand, such mixing have 



not been extensively studied further in any other model independent way. In this respect, it is 
important to derive model independent constraints in QCD that can be applied to constrain 
the phenomenological parameters in any model calculation. Here, we will construct the 
QCD sum rule for the vector-scalar correlation function and derive a constraint equation for 
their mixing. Such procedure have been previously used to constrain the model parameters 
in the momentum dependent part of the light vector meson self energy at nuclear matter 18 . 

One caveat in working with the scalar meson is the heated discussion on its dominant 
quark content, which could be a tetraquark[l9|, [lo|. In the lattice approaches, most of the 



round state mass of 
while some 
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calculations using a two quark interpolating field seem to predict t he g 
scalar particle in the isospin 1 channel to be above 1.3 Gev[2ll. I22I . |23| . 
predict it to be around 1 GeV|26_|, [27]. In contrast, most lattice calculations based on a 
four quark interpolating field current consistently predict the mass to be around 1 GeV for 



Some calculations based on two quark 
30 . So while the lattice calculations 



the f , 00 0, Q and around 600 MeV for the a [2 
current also give some consistent prediction for a [2 
seems to favor a tetraquark picture for the scalar, a two quark component is not ruled out. 
This is in a sense not surprising as the true scalar should have both a two quark and a four 
quark component, and the real question is what their relative contribution to the total wave 
function should be. 

In the QCD sum rule approach for the scalar meson, a previous work based on a two 
quark interpolating field[3l|, is known to be controversial [32]. This is so because the single 
instanton configuration is expected to contribute to the correlation functions between the 
scalar currents or between the pseudo scalar currents, and spoil the respective convergence 
of the operator product expansion (OPE). In the present analysis, we will investigate the 
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correlation function between the scalar and the vector current, where both interpolating 
fields are composed of a quark and an anti-quark. Thus, the OPE are free from the single 
instanton contribution and a reliable QCD sum rule analysis will be possible. As we will see, 
by investigating the scalar through a two quark current, and modeling the phenomenological 
side through a particle hole contribution with known couplings, we find that for both ao(980) 
and a meson, its overlap to the two quark current is less than 20% of the values typically 
expected from the usual meson. Hence, our result suggests that the scalar nonet have a 



dominant tetraquark component |19j|. Furthermore, our analysis confirms that the scalar 
vector mixing indeed takes place in the physical time-like region. While the a meson is too 
wide to be observable in the u meson emanating from nuclear matter, the will be in the 

P- 



II. OPE 



We start from the correlation function between the scalar and vector meson. 

n M ( 9 ) = i J d 4 xe^(T[J(x), J M (0)]) (1) 

where J a > a ° = \{uu± did) and J^< p = ^vq^u ± dj^d). 

In the vacuum, the correlation function is zero as there can not be any coupling between 
the scalar and vector current. However, in nuclear medium Eq.flT]) does not vanish, as the 
medium provides the necessary four momentum. 

The OPE of Eq.Q at nuclear matter can be obtained in the standard way [1, [H, 34]- To 



leading order in coupling, the first non- vanishing operator occurs at dimension 6. Taking 
the quark operators only, it is given as 



n „ = W jj_g /g __ as j ((uu)o{uru) n . m . + (dd) {dYd)n. m ), (2) 

for both the isospin 1 and channel. Here, (-) ((-) n .m.) denotes the vacuum (nuclear matter) 
expectation value. We have assumed vacuum saturation to extract the vacuum expectation 
value of the quark condensates. For a symmetric nuclear matter at rest, the expectation 
value becomes to leading order in nuclear density, 

((uu) (u^u) + (dd) (d^d)^j - (qq)oPN^°, (3) 

where pN,mN are the symmetric nuclear matter density and nucleon mass respectively. Here 
we have made use of the linear density approximation (-) n . m = (-)o + ^-(N\ ■ \N). 



III. PHENOMENOLOGICAL SIDE 

Consider taking the nuclear matter expectation value of the correlation function in Eq.(0Q). 
Since this correlator vanishes in the vacuum, it is just the nucleon expectation value times 
the density factor in the linear density approximation. If we saturate the intermediate states 
by hadronic states, the contribution from the ground states becomes as follows 

= ^A^mff - ?V7/)^hr. (4) 

2m,N q — m z s T V / q 2 — m v 
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Here, 

(0\J\S) = fs 

(o\r\v) = ffr, (5) 

where S, V is the scalar and vector meson, and e M the polarization tensor of the vector 
meson. Mg +N _> v+N is the forward scattering matrix element. 
Let us assume the following phenomenological Lagrangian, 

£ = [g v N lll T a N - -^-Na^Nd^ V* + g s NSr a N, (6) 

where we define r a = a a (a a : Pauli matrices) for the isospin 1 channel and r a = 1 for the 
isospin channel. Using this, the forward scattering matrix element has the following form. 

M" s+N ^ v+N = g v g s (g"» - q'q'h 2 ) g2 -^^ 2m N Nl„^fN 

- 9V/ 4^4 ( f.,)v/ m ^ (r ° )27V (7) 

Substituting Eq.([7j) into Eq.(jlj), one finds that Eq.(J3]) can be written in the following form, 

= n(g 2 , p ■ qf h ™ x [g^ - q»q»/ q 2 ) N^(r a ) 2 N. (8) 

For the nucleon at rest, the tensor part combined with the nucleon expectation value is 
proportional to \q\ 2 . This means that the vector scalar coupling vanishes when q —>■ 0. 

The sum rule constraint is obtained from identifying the W hen in Eq.fjHJ) to the corre- 
sponding OPE in Eq.(j2J). In the limit where q = (u, 0,0,0), we find the following sum 
rule. 

( 1607r V \ 1 



where the dots represent contributions from excited meson states. 

The first term in the right hand side of Eq.([9]) represents the contributions from the vector 
coupling in Eq.pl), while the second that from the tensor coupling. It should be noted that 
since we have calculated only the leading term of the OPE, the QCD sum rule will constrain 
only one parameter in the phenomenological side. 

IV. SUM RULES 
A. to — a mixing 

Let us start from considering the isospin channel. Here, the mixing of the uj will 
be to the a. As discussed before, many phenomenological and theoretical considerations 



■ = -\fsfvgs9v 
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lead us to believe that the dominant quark content of the scalar nonet is a tetraquark [It 
The previous QCD sum rule calculation for the isospin and 1 scalar meson with a two 
quark interpolating; current, seemed to reproduce the masses of a , fo to be of 1 GeV and 
degenerate 3~D] . However, direct instanton contribution was found to break the degeneracy, 
such that in the isospin 1 channel, the sum rule couple dominantly to do (1450), and in the 
isospin channel to a state with a mass smaller than 1 GeV, which could be the er(600)[35]. 
On the other hand, recent calculation with tetraquark interpolating currents do seem to 
reproduce the masses _of the scalar nonet correctly, if direct instanton contributions are 
taken into account I3j 
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The culprit for the complication in the scalar channel is the direct contributions from 
a single instanton j32|. In the present work, since we are considering the non diagonal 
correlation function between the scalar and the vector current, no direct instanton will 
contribute to the zero-modes. Moreover, since the ground sate will dominate the sum rules, 
we can get direct information on fg, which is the overlap of the ground state scalar meson 
to the two quark interpolating field. 

The mixing in the isospin channel is slightly complicated because the a width is very 
large and the u will also mix strongly with /q(980). On the other hand, we have calculated 
only the leading term in the OPE, and consequently can constrain only one parameter. 
Therefore, as a first approximation, we will take the o width to be small, neglect the con- 
tribution from /q(980), and assume that the u nucleon coupling is dominated by the vector 
part such that k = 0, as is phenomenologically motivated. We then perform the Borel 
transformation, which will suppress the contribution from / (980), and then compare the 
OPE to the phenomenological side. 

After the Borel transformation, the sum rule for the uo — a mixing becomes, 
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(10) 



The sum rule in Eq.fJTQl) can be used to constrain the whole coefficient f a f^g a gu}- However, 
phenomenologically g<7, fw ar e all rather well known. Therefore, we will use the sum rule 
to constrain the parameter f a . For the parameters, we take g w = gNNu = 11.5, g a = gNNa = 
9.4[!lj], m a = 0.55 GeV, and f u = 0.12 GeV 2 44]. The same value of can be obtained from 
the QCD sum rule calculation for the vector mesons|33] from which we also take f u = f p . 
For the parameters appearing in the OPE, we will take the following values, 



(qq) M=lGeV 

a(M 2 ) 



-(0.23GeV) 3 , 
Air 



91n(M 2 /(0.2GeV) i 



(11) 



Fig 1 shows the sum rule result for f a normalized by !(0|gry 5 r°g|7r ) = 0.13GeV 2 calcu- 
lated in the soft pion limit. One notes the sum rule has a plateau at Borel mass of around 
1.0 — 2.0 GeV 2 , at which the overlap is about 10% of that of the pion. This value fits very 
well to the picture that the scalar nonet is mainly a tetraquark and has a very small quark 
anti-quark component. Because the parameters for a have some uncertainties, we check the 
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sensitivity of our result to m a and g a . For m a uncertainty, we reanalyze the Borel plot for 
m a ~ 0.4 — 0.8GeV, and confirm that the change of the result is less than 10%. This seems 
surprising at first when looking at the individual terms in the right hand side of Eq. lflOl) . 
But it just reflects the fact that for the overall sum, the dependence on m a is weak, as 
can be seen from expanding the right hand side of Eq. fllOl) in 1/M 2 and noting that the 
leading term has no dependence in m a . For the uncertainty of g a , we consider the larger 
value of g a = 14 — 17 which is recently obtained 0,13. However, using these values leads 



to reducing our result of f a by about a factor of 1/2, which indicates the quark anti-quark 
component is further suppressed. 

Finally, we comment that the overlap f a to that of the pion overlap constant (= 0.13GeV 2 ) 
is not special as the corresponding overlap value for the vector meson is f p — — 0.12 GeV 2 
and f ai ~ 0.17 GeV 2 (3l| for the axial vector meson, and therefore, f a is much smaller than 
any of these values. 



B. p — a-o mixing 



The sum rule for the isospin 1 channel can be obtained similarly. In this case, however, 
the p coupling to the nucleon is dominated by the large tensor coupling and we will assume 



9p — 9nn p = 0. For the other couplings, we will take k = 14.2 and g 



(Hi 



9NNa 



2.8 



from the Bonn potential 41J. For the overlap constant for the vector meson, we again take 



f p = 0.12GeV 2 44j. The Borel sum rule for this case then becomes, 
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(12) 



The sum rule for f ao can be obtained by dividing the right hand side of Eq.(fl2l) by its left 
side apart from f ao . Fig. [2] shows the the value for f ao , normalized by the corresponding pion 
value. Unfortunately, there is no stable Borel region from which we can reliably determine 
the value. Adding the vector coupling of the p will not change much. The problem in 
this case, could be due to the nontrivial contribution from ao(1450), which is expected to 
be a dominant quark anti-quark state. In fact, previous QCD sum rule calculation with 
two quark current in the isospin 1 channel find the ground state mass to be around 1400 
MeV[36], suggesting that the a (1450) dominantly contributes to the sum rule. Therefore, 
in the sum rule in Eq. fll2l) . while the Borel transformation suppresses the contribution from 
the a (1450), its contribution might not be suppressed due to the large overlap / ao (i45o)- 
Unfortunately, the relevant couplings for ao(1450) are not known, and can not be subtracted 
out from the sum rule to obtain f ao . 

Another reason why the sum rule in Eq. (1121) is less reliable than that in Eq. fTlOl . can 
be seen from the phenomenological side. Apart from the couplings, the right hand side of 
Eq. ffl2l can be obtained by taking M A d/dM 2 of the right hand side of Eq. ffl0|) . Such deriva- 
tives tend to enhance the contributions from higher energy states, as additional factor of 
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resonance mass is multiplied in front of the exponential suppression factors. Therefore, while 
the contribution proportional to e _m<T//M is larger than that proportional to e~ mu ' M and 
that to negligible e -4 " 1 ^^ 2 in Eq. ffTOl) . the contribution proportional to e~ nia o^ M becomes 
less dominant in Eq. (fT2l . suggesting the importance of contributions from excited states. 

However, since we have calculated the OPE to the leading power corrections only, includ- 
ing additional terms in the phenomenological side becomes meaningless. Moreover, we will 
take the asymptotic value at higher Borel masses, where the approximation of taking the 
leading power correction becomes more reliable. At larger Borel mass, the curve in Fig. §2§ 
approach an asymptotic value for f ao at about 20 % of the corresponding pion value. This 
is then quite consistent with the case for the a case, and with the fact that the dominant 
quark content of the scalar nonet is a tetraquark. 

With all our result, we have shown that a consistent picture emerges from the QCD 
sum rule analysis where there is a nontrivial coupling between the ao and the p in the 
nuclear medium, whose strength can be reliably estimated with previously determined cou- 
pling constants. Therefore, such corrections should always be included in estimating the a 



contribution to the dilepton spectrum from heavy ion collision [16l. Il7 



V. SUMMARY 

We have derived the QCD sum rules for the vector scalar mixing at nuclear matter in 
both the isospin triplet and singlet channel. Since the phenomenological parameters are well 
known except for the overlap of the scalar interpolating field to the corresponding ground 
state scalar mesons, the sum rule can be used to calculate the value of the overlap. We find 
that the overlap in both the 1=0 and 1 channels are very similar but less than 20% of the 
corresponding value in the pseudo scalar channel. This result confirms that both the a and 
ao (980) and probably the remaining scalar nonets have a small quark anti-quark component 
and thus suggests a large tetraquark component. If the coupling of the scalar nonet to the 



nucleon is very large as in the estimates in [43, 1451 ] . the overlap would be even smaller. 

The mixing between different spin states originates from the fact that the nuclear medium 
provides a 4 momentum for the different spin states to couple. If the nuclear medium is not 
isospin symmetric, then there could be mixing between different isospin states also. So for 
example, in neutron rich matter, there could be coupling between different isospin and spin 
states, such as between the ui and the a (980). 

Experimentally, these mixing poses new challenges, as the extra peaks from mesons with 
different spins could be observable from the vector mesons emanating from the nuclear 
medium. As we have shown within the QCD sum rule analysis, a consistent picture of mixing 
between scalar and vector meson emerges, whose strength can be consistently estimated with 
previously determined phenomenological couplings. Therefore, while the a peak is too wide 
to be observable from the uo emanating from the nuclear medium, the cto peak will appear 
nontrivially in the p meson channel. 
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FIG. 1: Sum rule for f a = ^(0\uu + dd\a) normalized to ^(0|m7 5 u — di^d\ir°) . 
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Tables 



Isospin 


fv 


9S 


av 


K 


i j-. . c f s 7 . r „ „. from present work 

tj(U \wy°u— di~f J d 7r u ) 





0.12 GeV 2 


9.4 


11.5 





0.10 


1 


0.12 GeV 2 


2.8 





14.2 


0.2 



TABLE I: Parameters for the ui — a (p — ao(980)) mixing in the isospin (1) channel. 
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